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In this first paper we present a Lagrangian framework for the description of structure formation 
in general relativity, restricting attention to irrotational dust matter. As an application we present 
a self-contained derivation of a general-relativistic analogue of Zel'dovich's approximation for the 
description of structure formation in cosmology, and compare it with previous suggestions in the 
literature. This approximation is then investigated: paraphrasing the derivation in the Newtonian 
framework we provide general-relativistic analogues of the basic system of equations for a single 
dynamical field variable and recall the first-order perturbation solution of these equations. We then 
define a general-relativistic analogue of Zel'dovich's approximation and investigate its implications 
by functionally evaluating relevant variables, and we address the singularity problem. We so ob- 
tain a possibly powerful model that, although constructed through extrapolation of a perturbative 
solution, can be used to put into practice nonperturbatively, e.g. problems of structure formation, 
backreaction problems, nonlinear properties of gravitational radiation, and light -propagation in re- 
alistic inhomogeneous universe models. With this model we also provide the key-building blocks 
for initializing a fully relativistic numerical simulation. 
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I. INTRODUCTION 

General-relativistic analogues of the celebrated 
"Zel'dovich approximation" [74, 76, 79-81] for the de- 
scription of structure formation in the mildly nonlinear 
regime have been suggested previously, first by Kasai in 
1995 ([51]; for extensions to second-order perturbation 
solutions see [72], [73]). We shall put these works into 
perspective, as well as those by Matarrese and coworkers 
([61, 62]; for first- and higher-order perturbation 
solutions see [63], [64]), who discussed the relativistic 
analogues of the Newtonian equations in Lagrangian 
form. Croudace et al. and Salopek et al. discussed the 
Zel'dovich approximation in relation to spatial gradient 
expansion [39, 75], and Ellis & Tsagas [48], proposed a 
covariant form for the peculiar-motion corresponding 
to Zel'dovich's ansatz. The reader may also consult 
the seminal papers [44, 60] and [49] that address the 
application of the orthonormal frame approach to 
relativistic cosmology. 

In this series of papers we reinforce the Lagrangian 
point of view with full rigor within the framework of Ein- 
stein's equations, keeping the formalism as close as possi- 
ble to the Newtonian framework. In the present work we 
so obtain (i) a natural analogue (in form and in spirit) of 
Zel'dovich's model generalizing the approximation sug- 
gested by Kasai (loc.cit.), e.g., we obtain a quadratic 
form for the metric, useful for a realistic study of the light 
cone structure, together with nontrivial projected curva- 
ture and Weyl tensor approximations, including a nonlin- 



ear gravitational radiation part; (ii) general-relativistic 
Lagrangian equations that feature the Lagrange-Newton 
system of equations as a clear-cut geometrical limit (pro- 
viding an alternative to the set of equations derived by 
Matarrese & Terranova (loccit.) by insisting on a single 
dynamical field variable; and (iii) a covariant description 
of relevant kinematical and dynamical variables in the 
spirit of Ellis & Tsagas (loc.cit.). We shall also provide 
a number of useful details related to the basic equations 
and in particular to the electric and magnetic parts of the 
projected Weyl tensor that will be needed in forthcoming 
work. 

Paraphrasing a previous Newtonian investigation [14- 
16, 18, 42], we look for general-relativistic analogues of 
(i) the Lagrangian deformation gradient of fluid elements, 
(ii) equations that feature this variable as the only dy- 
namical one, and (iii) the corresponding first-order so- 
lution for perturbations at a FLRW background cosmol- 
ogy, which is then extrapolated into the mildly nonlinear 
regime according to a definition that we shall provide 
here. We shall restrict our investigations to the matter 
model "irrotational dust" . We obtain clear-cut answers 
to all of the above-mentioned points and discover a very 
close analogy between relativistic and Newtonian equa- 
tions and models. This allows us to easily transfer "New- 
tonian knowledge" to the relativistic stage. The success 
of the corresponding Newtonian approximation suggests 
that the relativistic version of Zel'dovich's model pre- 
sented here is a promising and possibly powerful one. 

This paper also aims at furnishing the basis for stud- 



ies of nonlinear perturbativc and nonpcrturbativc generic 
models that complement standard perturbativc studies 
and studies of exact solutions with high symmetry. This 
provides not only the basis for applications to structure 
formation in relativistic cosmology; future work will com- 
bine this approximation with exact evolution equations 
for the spatially averaged variables, [19, 20], yielding 
nonpcrturbativc models that are capable of addressing, 
e.g., the "backreaction problem" in relativistic cosmol- 
ogy, also as a possible source of "Dark Energy" , see the 
reviews [24, 25, 33, 38, 45, 54, 69] and references therein. 
Furthermore, these relativistic models open the door to 
other applications like, e.g., the understanding of nonlin- 
ear features of gravitational radiation, as well as of light- 
propagation and distance measurements in realistic inho- 
mogcneous universe models that all cannot be addressed 
within the Newtonian framework. Finally, since Newto- 
nian simulations are often initialized with the Zel'dovich 
approximation, future fully relativistic simulations can 
be initialized by the relativistic form of this approxima- 
tion as a first step. The ingredients needed for realizing 
the initial conditions architecture of such simulations arc 
provided here. 

A major motivation of this line of works is to put our- 
selves into the position to master a possible paradigm 
change in cosmology that entails the need for inhomo- 
gcncous relativistic models. Curvature effects may play 
a key role in accessing the interpretation and high- 
precision determination of cosmological parameters in 
the near future. For example, the averaged spatial 
scalar curvature may evolve differently from a constant- 
curvature homogenous model, starting with a small cur- 
vature as furnished by cosmic microwave background ob- 
servations and producing an effective negative curvature 
in the Late Universe [24, 27, 33, 68], bringing geomet- 
rical and topological features into the fore. Another is- 
sue is the interpretation of cosmological parameters that 
may be affected by curvature and Riemannian volume 
effects when comparing averaged variables in an inho- 
mogeneous geometry with averages on a Friedmannian 
template space (e.g., [26, 57]). Finally, note that New- 
tonian cosmologies require periodic boundary conditions 
for any model of structure formation [30] , which can be 
relaxed in a relativistic setting (see e.g., [21, 24, 25]). 

We proceed as follows. In Section II we recall the 
Newtonian derivation of the Lagrangian equations to- 
gether with Zel'dovich's approximation in terms of a 
first-order Lagrangian perturbation solution. In Sec- 
tion III we paraphrase the Newtonian derivation within 
general relativity, give a compact analogous formulation 
of Einstein's equations using Cartan's coframes, and dis- 
cuss first-order perturbation solutions. Section IV de- 
fines the general-relativistic analogue of Zel'dovich's ap- 
proximation and discusses it in full detail by function- 
ally evaluating relevant variables including geometrical 
fields. Section V proposes tests of the extrapolation into 
the nonlinear regime, discusses relations to the singular- 
ity problem, and highlights the main findings including 



follow-up prospects. Appendixes are dedicated to alter- 
native formulations of the governing equations and an 
example for the proposed approximation. 



II. LAGRANGIAN THEORY OF STRUCTURE 
FORMATION IN NEWTONIAN COSMOLOGY 

In this section we briefly recall the logical structure of 
a derivation of Zel'dovich's model within Newtonian cos- 
mology. Thereafter we shall contemplate on Zel'dovich's 
original suggestion and his extrapolation idea in order to 
prepare ourselves for the relativistic setup. 



A. The Lagrange— Newton— System 

In the framework of Newtonian gravitation the field 
and evolution equations governing the motion of self- 
gravitating dust form a closed system in the Eulerian 
picture, consisting of the Eulerian evolution equations 



dtv = - {v ■ \/) V + g , 


(1) 


dtQ = -V • [qv) , 


(2) 


and the linear gravitational field equations 




V X .g = d. 


(3) 


V -g =: A - AtiGq . 


(4) 



We call this system of equations the Euler-Newton- 
System (ENS). Here, as usual, g is the dust's density, 
G the gravitational and A the cosmological constant. 

Now we perform the transition from the Eulerian to the 
Lagrangian picture. Then, as we shall see, the trajectory 
field X = f{X,t) - defining the coordinate transforma- 
tion at a fixed time, or a time-dependent diffeomorphism 
- will be the only dynamical field variable that remains in 
the transformed equations, where X'^ are the Lagrangian 
coordinates, comoving with the fiuid, that are defined as 
to coincide with the Eulerian ones at some initial instant 
of time. The field f{X,t) measures the deviation of a 
fluid element's position at some time t from its initial 
one, and its Lagrangian gradient (/'i,) measures the vol- 
ume deformation of fluid elements, where a vertical slash 
is used to denote partial derivative with respect to La- 
grangian coordinates. Upon introducing the trajectory 
field f{X,t) we implicitly solve the Eulerian evolution 
equations by 



f ; 9^ f 



Q 



J>0 



(5) 



with the initial density field giX), and J the Jacobian 
determinant of the transformation from Eulerian to La- 
grangian coordinates (using the functional determinant 
notation in the first expression). 



T _ '^l/ 1 J 1 J ) _|_ Imn fi fj fk 

~ d{X^,X^,X^) ~ 6 ■''' '' 1™ I" 



(6) 



The Eulcrian field equations then assume the form of La- 
grangian evolution equations, if the field strength is ex- 
pressed through the trajectory field as above. The result- 
ing system of equations, the Lagrange-Newton-System 
(LNS), takes the following form ([31] for A = and [14] 
for A^O): 



J{f\l\f)+cyd. 



0, 
A J 



AtiG'q , 



(7) 
(8) 



where J^ denotes the functional determinant of the ex- 
pressions in brackets. Other forms of the Lagrange- 
Newton-System may be found in the review [42], and 
in Appendix A. 

Finally, we introduce the Newtonian tidal tensor £ij , 

= ^j^,kiJ{rj\f) + ^(47rG4 - A)<5^ , (9) 



2J 



J 



where we have inserted the field equation (4) in the sec- 
ond line. In terms of this form of the tidal tensor, written 
as a functional of /, we can express the LNS through the 
symmetry conditions on the tidal tensor £ij : 



^M = 



and 



E\ = 



(10) 



furnishing the four Lagrangian evolution equations for 
the three components of the trajectory field. 



B. Derivation of a first-order scheme and 
Zel'dovich's approximation 

Now we proceed by linearizing Equations (7) and (8) 
at a reference background with respect to the deviations 
from this background: we assume the only variable / to 
be a superposition of a homogeneous and isotropic back- 
ground deformation fH{X,t) = a{t)X and an inhomoge- 
neous deformation field p{X,t), i.e. 



fiX,t)^ait)X+piX,t), 



(11) 



where for convenience a{to) := 1 and p{X,to) = 0. It is 
sometimes useful to introduce the scaled trajectory field 
(f = F(X,t) = f{X,t)/a{t) to describe motions in a co- 
ordinate frame q that is comoving with the background 
solution. Correspondingly, we may also introduce the 
scaled deviation field P{X, t) = p{X, t)/a{t) that will be 
useful when comparing with the relativistic setting. 

A homogencous-isotropic deformation separately 
solves the LNS. This yields Friedmann's expansion law 
as a first integral, with background density qh = QHCi~^ , 



H'- 



SttGqh + A const 



(12) 



The first -order system of equations, which has to be 
solved, is 



a Vo X p— aaVo x p = , 



(13) 



a^Vo • l>+ {2da - a^k) Vq • p = -47rG [g - qh) ■ (14) 

Here Vq denotes derivative with respect to the La- 
grangian coordinates. If we insert the field equations 
(3) and (4) at initial time, we are able to express the 
source term in (14) by the divergence of the initial field- 
strength perturbation p(io). Now, we split the perturba- 
tion field p into a longitudinal part p^ and a transverse 
part p"^. The resulting equations are for the transverse 
(divergence-free) part: 



r --r =0, 

a 
whereas the longitudinal (curl-free) part obeys: 

1 -v, 



p"+ 2--A p" = -p-(to) 



(15) 



(16) 



Using the deviation field P{X,t), Equations (15) and 
(16) take the form 



P^ + 2i/P^ = 0, 



P^ + 2HP^ - AttGqhP 



(17) 
-W{X) , (18) 



where W{X) = P^{X, to)+2H{to)P^{X, to) is the initial 
peculiar-acceleration field. A more detailed derivation 
together with a general solution to these equations can 
be found in [15]. We wish to refer the reader also to 
[42] where the general equation and solution schemes for 
perturbations at any order can be found. 

As a special case of the general first-order solution we 
obtain the "Zel'dovich-approximation" , when we restrict 
the peculiar-velocity and peculiar-acceleration fields at 
some initial time io by the "slaving condition" : 



u{X,t) = w{X,t)t ; t = to 



where 



u = aP , and w ^ u + Hu ~ 2aP + aP 



(19) 



(20) 



The condition (19) is then preserved in time. 

In the case of a spatially flat background without 
a cosmological constant the "Zel'dovich-approximation" 
reads 



NZA 



F{X,t)=X + - 



^'- 



P^(X,to)to. (21) 



For general backgrounds including a constant-curvature 
term and a cosmological constant, sec [12]. 

As we shall discuss in more detail later, Zel'dovich sug- 
gested to extrapolate this trajectory field into the mildly 



nonlinear regime, so that the nonhnearly evolved density 
can be calculated through its exact integral 



NZA, 



QHJt) 



giX)/''^^JF{X,t), 



(22) 



where Jp = det(i^^ ) is evaluated for the comoving tra- 
jectory field (2f ). His motivation was that this expression 
for the density, if linearized at the background, coincides 
with the linearized solution for the density in comoving 
Eulerian coordinates q, while the nonlinear expression is 
capable of describing a continuum that develops caustics 
in a finite time, in form similar to the rectilinear motion 
of an inertial continuum [82] (for further discussions of 
this extrapolation idea and its subsequent developments 
see [14]). 

We summarize the logical structure of the derivation 
of Zel'dovich's approximation: (i) the basic system of 
equations furnishes a closed system for a single dynami- 
cal variable, the deformation field /, or the deformation 
gradient /V ; (ii) introducing a split into a background 
deformation (the Hubble fiow) and a deviation field, we 
exploited the fact that we have only to linearize in the de- 
formation field and not e.g. in the density deviations as in 
the Eulerian picture; (iii) Eulerian fields, e.g. the density 
field, but also others, can then be evaluated as functionals 
of the linearized perturbation and so provide nonlinear 
expressions as an extrapolation into the mildly nonlinear 
regime (i.e. up to shell-crossing singularities develop, af- 
ter which the transformation of the Lagrangian function- 
als back to Eulerian space is no longer regular). The fur- 
ther restriction of initial data is not mandatory so that, 
in principle, we can use this extrapolation idea also for 
the general first-order solution including vorticity. The 
functional for the vorticity is given by Cauchy's exact 
integral [15]: 



a^Jp 



nEEOj{X,to) 



(23) 



ric as a quadratic form of the deformation field will re- 
main a quadratic form in this approximation. We are so 
able to keep highly nonlinear information encoded in the 
functional dependence on the perturbation variable (e.g. 
the exact density integral, the Ricci and Weyl curvatures 
etc.), while their solution is explicitly expressible in terms 
of constraint initial data and known time-dependent co- 
efficients. While Zcl'dovich and his coworkers mainly ex- 
ploited the nonlinear functional dependence on the de- 
formation in the density field, we here wish to apply this 
logic to all functionals of interest. As emphasized pre- 
viously, this strategy is only applicable if the governing 
equations form a closed system for the deformation vari- 
able alone. 



III. LAGRANGIAN THEORY OF STRUCTURE 
FORMATION IN RELATIVISTIC COSMOLOGY 

In this section we shall introduce the coframe field be- 
ing the generalization of the Lagrangian deformation gra- 
dient of Newtonian cosmology. In the general-relativistic 
case the deformation of fiuid elements is no longer inte- 
grable, i.e. instead of the basis dx°' = f^.AX"^ we have 

to consider a non-exact basis r/" = rj'^^dX^ . While the 
linearly transformed (Lagrangian or local) basis in the 
Newtonian case derives from three functions (the com- 
ponents of the trajectory field), here the linearly trans- 
formed local basis (here viewed in the cotangent space 
at a point of the manifold) involves nine functions (the 
coefficients of the set of coframe fields); hence we have 
to find at least nine evolution equations. (We use latin 
letters a, &,••• = 1, 2, 3 as counters in order to distingish 
them from coordinate indices i,j, • • • = 1,2,3; through- 
out the paper || denotes covariant derivative with respect 
to the 3-metric with a symmetric connection, whereas | 
denotes partial derivative as before.) As in the previous 
section the three spatial deformation one-forms will be 
the only dynamical variables in our setup. 



The strategy to find the corresponding 
relativistic approximation 



A. The Lagrange— Einstein— System 



According to what has been said above, a general- 
relativistic analogue of Zel'dovich's approximation has 
to aim at (i) writing Einstein's equations in terms of a 
system of evolution equations that all feature a single dy- 
namical field variable corresponding to the Lagrangian 
deformation gradient; (ii) reducing constraint equations 
to constraints on initial data where possible; (iii) find- 
ing the general first-order solution of the system of evo- 
lution equations for the deformation variable, and then 
(iv) employing Zel'dovich's extrapolation idea to func- 
tionally express other variables in terms of the single 
perturbed deformation. It is clear that such a strategy 
results in a nonperturbative approximation of relevant 
field variables. For example, the resulting spatial met- 



Restricting the matter model to "irrotational dust", 
the simplest spacetime foliation is given by a family 
of flow-orthogonal hypersurfaces with induced 3— metric 
coefficients gij in the comoving and synchronous metric 
form, 



(4) 



g 



dt ® dt - 



(3), 



(3) 



^ = g,,dX'®dX\ (24) 



where A"' are Gaussian normal (Lagrangian) coordinates 
that are constant along flow lines (here geodesies). (The 
proper time-derivative is equal to the coordinate time- 



derivative, u^TT-n- 



— , below written with an overdot; 



dt 



Greek indices label spacetime and Latin ones spatial co- 
ordinates.) In our foliation we define the coefficients of 



the extrinsic curvature as usual, 



1 



r^ij — ^9ij 



(25) 



In the following we replace this quantity by the expan- 
sion tensor coefScients 6(j, •> = ^^ij together with the 
symmetry condition 6rj,i = (required in our foliation). 
In this frame the Einstein equations take the well-known 
form which has been introduced by Arnowitt, Deser and 
Misner [10], consisting of the evolution equations 



6' 



g + Sg^O (26) 

+ OQ'j = {AnGg + A) S'^-R) , (27) 



where Rij are the Ricci tensor coefficients corresponding 
to the 3-metric, and g the dust density field, completed 
by the constraint equations 



e,, - e, 







i? + e^ - e'^^e' fe = mnCg + 2A 



e' 



J life 



l^ k 



0, 



(28) 
(29) 
(30) 



The latter equations are the Hamilton constraint and the 
momentum constraints. 

In the (3 + l)-split the raising and lowering of indices 
docs not commute with the time derivative, i.e. for any 



tensor Tij we find 



t,^i9^kT)y^g.kT) + 2e,,T'^.^ 



(31) 



In the same way the covariant derivative and time- 
derivative do not commute. To handle this we derive 
the following useful relation for the symmetric connec- 
tion coefficients (Christoffel symbols): 



-L fe/ — *-^ k I 



^\U - 9''Q 



mi 



(32) 



Using it we are able to rewrite the time-derivative of the 
3-Ricci tensor coefficients through spatial derivatives of 
the expansion tensor coefficients, 



R, 



e' 



«lb1|fe 



e' 



iMk 



-e 



life 
u II fc 



-e 



\i\\j 



(33) 



We can thus recast the evolution equation for the expan- 
sion tensor (27) into a form that only features expressions 
built from the expansion tensor and its derivatives (apart 
from the exactly integrable source g): 



e' 



ee' 



e 



i II fc 






2e\.e^ 



2ee'fce'=j- 



e 



k II J 



e* 



II fc 



J life ' ^ \\j ^ 1 life ^ feiij 

2 (47rG£- -h A) e*j- - 47rGee<5'j . 



(34) 



We note already that the trace of the expansion tensor is 
G = J I J ., with J given below, Eq. (38), so that we can 
immediately solve the continuity equation g-|- Gg = by 
integration and get the general integral in analogy to the 
Newtonian case: 



It is also useful to note that, in view of the momentum 
constraints (30), we have g^^R^i = and therefore 



R 



-m\R^ 



(36) 



1. Using Cartan's coframe fields 

Introducing Cartan's coframcs rj" which define, up to 
rotations, a noncoordinate basis of 3-dimensional space 
we rewrite the spatial part of the metric as 



(3) 



g = SabV" ® V'' 



gij = Sabif^rf^j 



(37) 



Noncoordinate indices are raised and lowered by 5^^^. 
Our choice is to simplify calculations by putting all the 
information on the initial data into the coframes by 
rj°'i{to) = 'n"'i- However, one could choose the more gen- 
eral orthogonal (and not orthonormal) matrix G^jj, in- 
stead of (5^^; then the coframe would take a simple form 
at some initial time, '^'r]°'i{to) = (5° , which would formally 
come closer to the Newtonian description. We shall keep 
the standard definition throughout this first paper, but 
we shall come back to the other choice in forthcoming 
papers. Chandrasekhar [34] discusses circumstances in 
which such a more general choice is useful. 

Throughout this paper we define the Lcvi-Civita- 
tensor density by ^i-^i^i^ ~ (—1)^, where P is the sign 
of the permutation (1,2,3) — >■ (^112^3) and e^ j j = if 
any two indices are the same. 

The determinant of the transformation between the 
coordinate and noncoordinate basis is given by 



J 






(38) 



The (inverse) orthonormal vector basis is described by 
the triads (frames) ea = e^d/dX"^, which can be ex- 
pressed in terms of the coframes as follows. 



1 



Ca = ^eabce''''v\v"i 



(39) 



Thus, the coefficients of the inverse metric take the form 

ij _ cab i j _ 

9 -0 P^a^b - 2 J2 



<^ce%e''-'e^""^^fc^'zry^™r;/„ . (40) 



We rewrite the vanishing of the covariant derivative of 
the metric g Ui^uu = using frames and coframes and 
get our first two constraints (for the second one we apply 



^ [»lbl -'' [W 



^\m)^ 



eaVfc||,=0, 



(41) 
(42) 



g=^ g 



1 



(35) 



Since J = y^ in our choice of coordinates, where g is the 
determinant of the 3-metric, we have to keep in mind 



that J is not a scalar but a scalar density, i.e. its covariant 
derivative vanishes according to (41), 



J ^ 2J 



= ^T^tabce 77 ,„||,,77 kV i^e^V k\\r = , (43) 



but its partial derivative with respect to the spatial co- 
ordinates does not, 



Ji 



7 = e.S\|, = r\, ^ . 



(44) 



It is also important to note that, contrary to the New- 
tonian case, the initial transformation determinant does 
not equal to one, J 7^ 1 (as it would do for S^i^ -^ G^f^). 
The invariant volume element then is described by the 
tensor 



ikl 



Je^ki and e'^^ = -e''=' 



(45) 



Since we use an orthonormal noncoordinate basis, we 
simply have e^^^ = e^j,^. We now turn to rewriting the 
ADM equations in terms of the coframes only. The ex- 
pansion tensor coefficients with one index inverted, 0* , 
will provide the closest analogy to the Newtonian case: 

©^ - ^,^: = ^^abce'''''h'',v\^\ , (46) 

and the symmetry condition (28) becomes 

Qm - S,,r,^[X,] - . (47) 

Now, the three latter constraint equations have become 
evolution equations. 

We define the 3-Riemann curvature tensor via the 
commutation relation of second covariant spatial deriva- 
tives, V^i\k\\i - V^i\i\ik ^ Rkirv^m, so we get: 

RkU3 -= 5ab{^l%\\k\\i - 'Allillfe)^^ ' 
and, finally, 

R'jki ^^a^'n°']\\k\\i ~^°j-||iiife) 

= ^^a.ce™" (^f.iifcii, - V^'m^) ^ ™^'" • (48) 
Contraction yields the Ricci tensor, i.e., Rij = R'' j^,-, 

i?., = ^^a.ce^'"" (77^11^11, - ^",||,|„) V\,V\ ■ (49) 

Note that the by simplifying this expression using (39) 
and the identity 



we get the alternative expression 

R - 5 hin" "'' -n" "'""In'' 

^ij — Oab(ll k ||j '/ fc||i )'l J 



(50) 



Finally, we find for the mixed 3-Ricci tensor coefficients, 
expressed solely with the help of coframes. 



R j -dab[V k -V k )V j 



(51) 



or, alternatively, with the help of frames (that often sim- 
plifies calculations). 



i?%. = <5"V(e,V||,-e,^||,||,). 



(52) 



Contracting the Ricci tensor, i.e., R = R^'k-, we obtain 
for the scalar curvature in terms of coframes. 



i? = <^afc(rA"'"'-r/\"''^"')rA, 



and, with (52), in terms of frames, 



ab ^ k(^ I 



i? = <5'^V(e 



; 



b \\l\\k ^6 ||A:|| 



l) 



(53) 



(54) 



Using the coframe as the single dynamical variable the 
ADM equations become: 



1 



^abc'^ 



'''r^v'kV'i) • = (47rGJ^ + A j) ,5'^. - Jff^- , (55) 



and the set of (former) constraint equations become 

SabV^v',] = ; (56) 

^abc^""''rrnV\v'i = ^Q^GJg + 2AJ - JR ; (57) 

{eabce"''r,vWl\, = i^abce'^'v^^vWlX, ■ (58) 



The first of these equations, (56), arises as the time- 
derivative of the symmetry condition for the expansion 
tensor, (57) comes from the Hamilton constraint, and 
(58) represent the momentum constraints. We still use 
the covariant derivative and the Ricci tensor in this 
(overdetermined) system of 13 evolution equations for 
the 9 components of the deformation coefficients 77"^, but 
do that only for the sake of readability. It is possible to 
express these in terms of the coframes only, too, as done 
above for the Ricci tensor and the scalar curvature. (We 
only have to make sure that all these equations including 
the covariant derivatives could be expressed in terms of 
the coframes only.) In this sense, an interesting form of 
(55) is the following: 
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V^^ + -ebcd^"''''v\v'',n]V%v''i - ( ^^G^ + A ) r;- 
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bcd'^ ["n"^ k\\i\\m ^ V"^ k\\m\\i)V I 



(59) 



Here the equation is solely expressed in terms of the 
coframes and their time and spatial (covariant) deriva- 
tives. 

We now condense this system of equations in a more 
compact form. First, we get another convenient form of 
(55) when we split it into its trace and trace-free parts. 
With the Hamilton constraint the trace takes the form 



eabce^"^°^\^A = AJ-4^GJe, 



(60) 



which is Raychaudhuri's equation. The trace-free part 
on the other side is 

2\^abc'^ V jV kV I o^afcc^ V mV kV l" jj 

-i-(c c^^'-'h'^ n'' rf- — —f f^nkl-a -b „c H \ 
+ \^abc'^ V jV kV I n^abc^ V mV kV l" jj 



= -J-) > 



(61) 



where r^ are the coefficients of the trace-free part of the 
Ricci tensor to be calculated from Eq. (50) and Eq. (53). 
We can now recast the set of equations (50), (50), (56) - 
(61) into what we call the Lagrange-Einstein System for 
dust (LES). 



2. Summary: the Lagrange-Einstein-System, 

The foUowing (over-determined) system of 13 evolu- 
tion equations for the 9 coframe coefficient functions is 
equivalent to the ADM set of equations for the matter 
model 'irrotational dust' (recall that this latter restric- 
tion implies that Eq. (62) already holds for the first time- 
derivative of the coframes): 



\^abc^''''v\v\ri\ = A J - AnGJ'Q 



(62) 
(63) 



eafcce"''^ m'7\^'i = le^GJe + 2A J - JR (65) 

\{^abc^'^'v'',v\v\ - \^abce"'''v''mV\ri\5\) 
+ {^abc^''''v'',v\v\ - \^abce"'''v''mV\ri\5\) 



-Jt\ 



(66) 



This Lagrange-Einstein-system is a system of equations 
described solely in terms of the coframes. We did not ex- 
plicitly insert the trace-free part and trace of the 3-Ricci 
curvature into the above equations because the resulting 
equations are tedious to read. In principle it can be done 
with the equations given in this section. Of course, the 
covariant derivative can also be expressed in the coframes 
language by calculating the Christoffel symbols with (37) 
and (40). 



The Weyl tensor is defined as the trace-free part of the 
4-Riemann curvature tensor, 

C'\, - '''R'\, ^K'''< + 1<S%(''R- (67) 



It has 10 independent components and thus carries all 
the information of the system. The Weyl tensor satisfies 
all of the symmetry conditions of the 4-curvature tensor, 
and in addition is trace-free over any two indices. It can 
be irreducibly split into two parts, called the electric and 
magnetic parts. Both parts are symmetric, trace-free 
tensors and have five independent components each. 



E„ 



C, 



flKl/X 



and H,. 






(68) 

After the (3+ l)-split the electric (tidal) part E"^ ^ of the 
Weyl tensor and its magnetic part iJ' take the following 
forms: 



E^ 
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e',e^. 
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(70) 



We also note the useful expressions 



ie. 



k^ ] 

k 
J 



of^ k^ 1° 3 



(71) 



where aij are the components of the shear tensor. 

We rewrite the above parts of the Weyl tensor by fully 
expressing them through coframes, 

E', = -i^^abce-'^r^ri^v'^i ^(4^G^ - a)<5'^, (72) 



H\ 



2J 

1 

1 



J 



'5a5e'''(^n-||;'A+'7",'/\||0 



(73) 



We infer that the projected electric part of the Weyl ten- 
sor yields a direct generalization of the tidal formulation 
of the Lagrange-Ncwton-Systcm (10): 



^[.fc] = 



£;* 







(62) 
(63) 



(74) 
(75) 



3. Form,ulation with the Weyl tensor 

We are now going to reexpress the above Lagrange- 
Einstein-System in terms of parts of the projected Weyl 
tensor in order to furnish the analogy with the tidal for- 
mulation of the Lagrange-Newton-System (10). Here, 
the electric part of the Weyl tensor plays the role of the 
tidal tensor of Newtonian theory, whereas its magnetic 
part carries additional information and describes gravit- 
omagnctic effects. 



whereas the magnetic part reproduces the momentum 
constraints and again the symmetry condition for the 
time-derivative of the expansion tensor: 



H, 



ik] 



Tjk 
ri I. 



(64) 
(62) 



(76) 
(77) 



At this stage the symmetry conditions on the electric and 
magnetic parts of the Weyl tensor do not cover all of the 
equations of the LES: the electric part of the Weyl ten- 
sor fully covers the 'electric part' of the LES, Eqs. (62) 



and (63) - which we need for the translation of the New- 
tonian approximation -, while the magnetic part of the 
Wcyl tensor just captures part of the 'magnetic part' of 
the LES, namely Eqs. (64), i.e. in total 7 equations. 
We believe that another form of the magnetic part could 
eventually provide a symmetric formulation of the whole 
system, but wc did not succeed to find it. 



4-. A geometrical Newtonian limit for spatial deformations 

It is easy to confirm that wc obtain the LNS as the 
following geometrical (spatial) limit to the Lagrange- 
Einstein-System of equations. The transition from gen- 
eral coframe coefficients ?/°j to those of the integrable 
(Newtonian) form 



V 



^^" 



/l 



(78) 



below) has no nontrivial Newtonian counterpart, as it 
always vanishes in the geometrical limit defined above. 



Nrjl _ 1 r 



Amn fa ( hh 



f%h^^r 



= 0. 



(84) 

This geometrical limiting procedure is spatial and, there- 
fore, does not involve a limit c — > oo. The light cone 
structure is simply not seen within the 3— space by co- 
moving observers (the Lorentzian structure appears in 
time-direction only). Note that this limit explicitly 
demonstrates that a Minkowskian limit (consisting now 
of this geometrical limit and c — ^ oo) in the comoving- 
synchronous slicing of spacetime is well-defined. In 
a post-Newtonian formulation, the Minkowskian limit 
leads to the Eulerian form of the Newtonian equations, 
while in this setting it leads to their Lagrangian form (see 
also [41, 43]). 



directly transforms the equations (62), (63) into (7), (8). 
This is particularly easy to see using differential forms 
as done in Appendix A. We also find that, as expected, 
the electric part of the Weyl tensor reduces to (minus) 



the tidal tensor of the Newtonian picture, Ei^ 



-£ 



The spatial line element then takes the well-known New- 
tonian form 

(3)g -> \^dX^^dX^ ^d,jy)^dX''®dX^, (79) 

i.e. a Euclidean line element that was transformed us- 
ing the transformation /. The basis vectors (frames) or- 
thonormal to the coframes in the Newtonian limit are 
/i ^ = ^e^j* (where a comma denotes derivative with re- 
spect to Eulerian coordinates). They obey f^ijhl = Sj 
and take the following form. 



hl = 



ikl fh fc 
WnJ abc^ J \kJ \l 



(80) 



In this limit the connection coefficients reduce to the 
inertial force terms 



N 



^ kl — J \kl'^,a 
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imn ra 



2Nj-ab 



/|fe//|™/V^O- (81) 



Although the Christoffel symbols do not vanish because 
in the Lagrangian picture of the Newtonian equations we 
adopt curvilinear coordinates, we can easily verify that 
i?j =0 as required and 
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„a _ fa 

V k\\i — J I 
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^ klT \i 



0. 



(82) 



However, since covariant derivative and time-derivative 
do not commute, we have nonvanishing 

Vfe||, = ./%,-/1A^7V,- (83) 

Combining the latter equations wc immediately confirm 
that the magnetic part of the Weyl tensor (see Eq. (73) 



B. Derivation of a first— order scheme 

To derive a first-order perturbation scheme wc choose 
a fiat homogeneous and isotropic background with some 
initial perturbation thereof, P°(X, io) = P^'iiX) (hence- 
forth, we omit the vector symbol over X for notational 
ease). With this choice of coordinates, the coframe coef- 
ficients take the form 

T?'', = a{t) [r, + P^^iX, t)] , and r^ ^ -v". , (85) 

a 

where a{t) is the usual scale factor and P'\{X,t) the in- 
homogeneous deviation (perturbation) field with respect 
to which we shall linearize the equations (for notational 
ease we shall not write ^'^■'P'^j). We call 77° the "peculiar" - 
coframe. At some initial time we assume a(io) = 1 and 
P'llX) 7^ 0. The initial perturbation cannot be set to 
zero because that would flatten the initial metric and 
suppress any metric evolution, as pointed out by Matar- 
rese & Terranova [61] as well as Russ et al. [72]. 

Thus, the initial 3-Ricci tensor is not equal to zero. Re- 
member that, generally, for a homogeneous and isotropic 
background within a space of constant intrinsic scalar 
curvature —6k/a^ the zeroth order of (55) is ^P%- = 
—2-^S'^j, see for example [51], whereas we have ^R^j ~ 0. 
However, since we can choose appropriate initial pertur- 
bations (i?*j ^ 0!) to describe the space we want, this 
choice implies no restriction of generality. 

The perturbation P° only appears summed over the 
noncoordinate index in the equations, so we introduce 
the following tensor coefficients and their trace: 



J a 



P" and 



P^P\ = dJ^P\ , 



(86) 



and use this notation throughout the remaining part of 
this section to make reading more convenient. (Note that 
within the first-order scheme wc have two true tensor 
indices here.) 



Field variables and equations in the first-order scheme 



Up to the first order the spatial metric takes the form 



(l)...-.2 



ik 



P,. + P 



]i) 



(87) 



(Recall that we here aim at strict linearization, not to 
be confused with the relativistic form of Zel'dovich's ap- 
proximation that we shall define below). 

The determinant of the transformation from coordi- 
nates to the noncoordinate basis becomes 

(i)j = a^{l + P) and (^'j = f 1 + pW 1 . 
The first-order Christoffcl symbols are 



(i)p< 



kl 



P 



(k\l) 



P ' 
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P, 



\^^ 



{kl) 



and in particular, we find ^^'r\j = P:- . 

To begin with the first-order LES equations, let us 
have a look at the symmetry condition (56). Straightfor- 
ward calculation up to first order yields 



P, 



M 







P, 
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P, 



M 



(89) 



We now derive the first-order expressions for the co- 
variant derivative of the coframe coefficients and its time- 
derivative, 



(1). 



a V k\\l 



a P 



lk\i] 



P, 
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{kl) 
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{k \i] 



Thus, we find an expression for the momentum con- 
straints (58) up to first order, which are 



pk 



■fe|j] 







P' 



[k\i] 



P' 



[m 



(90) 



Assuming the Hamilton constraint (57) holds for the ho- 
mogeneous background, that is, for vanishing perturba- 
tions, we get 



3ij2 = 8nGgH + A , 



(91) 



where qh = Qho/ o?' is the homogeneous background den- 
sity. Then the first-order Hamilton constraint reads 



HP + ^.ttGqhP 



iWp. 



(92) 



Next, we determine the Ricci tensor and its trace up to 
first order. Since the Christoffcl symbols (88) are of order 
©(P), their product is always of second or higher order, 
and we find 



(^)i?, 



(l)pfc __ {l)yk 

p _ p \^ 
^\v ^ij \k 

V(tk) \j + ^{jk) \i 



(93) 






2^K.) 



where we used the two constraint equations (89) and (90) 
above to express most of the terms by means of the initial 
perturbation field. 

Obviously, the first-order trace is solely dependent on 
the initial perturbation, 

so the first-order Hamilton constraint finally yields 



HP + AttGqhP = \P 



k \l 

[k \l] 



(95) 



Now, we address the six evolution equations of the LES 
(55). The homogeneous equations reduce to the trace 
and, with (91), yield Friedmann's acceleration law 



3- = -47rG£)H + A, 



(96) 



where qh ~ QHo/d^ again is the homogeneous back- 
ground density. The first-order equations then are 

P^ + ZHP'^ + HP5'^ + AtiGqhP&'j = -*^^i?^ , (97) 
and using the Hamilton constraint (92) we get 



p^ 



3HP\, - \hP5', - ^qhPS'j = -('^T^ , (98) 



3 



where '^^V*- = '^^^P* — i^^^'PJ' ■ represent the trace-free 
part of the 3-Ricci tensor. Hence, it is convenient to 
write these equations in their representation (60) and 
(61), i.e. split into trace and trace-free parts. The first- 
order Raychaudhuri equation takes the following form: 



P + 2PP - AttGqhP = . 



(99) 



We get another very convenient form of the scalar 
equation up to first order, when we take the trace of 
(97) and replace 3iJP by applying the Hamilton con- 
straint (i.e. instead of replacing the Ricci scalar as done 
in (98)). 

Since the first-order trace ^^■'P only depends on the ini- 
tial perturbation field, we subtract the resulting equation 
at t = fo from the general one to get rid of it. Hence, the 
alternative form of the trace equation is 



P + 3PP 



:{P{to)+3H{to)P{to)). (100) 



Note that this equation is equivalent to (99) (see the gen- 
eral solution in the next section). However, because of 
the following equations, this equation seems to be the 
more natural choice even if (99) is well-known from New- 
tonian theory. 

Now, the trace-free part of the set of evolution equa- 
tions (61) is 



H'^- + 3Pff ^. 



-Wt'- 



(101) 
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where H* ■ = P* — ^PS^ ,. This last set of equations gov- 
erns gravitational radiation, see also [61]. This is made 
clear if we express ^^V' ■ by (93) and (94) and - anal- 
ogous to what we did with the alternative form of the 
trace equation - choose the more convenient form 



n. 



3HU,, + 4 (n,, , - n, J^,) = 4c,-,- , (102) 



ij \kl 



t2 U ' 



where Cij are constants depending on the initial pertur- 
bations, 



c,, = ri,,-(io) + iH{to)n,Ato) + n|,,(to) 



n 



ifc 

ij I k 



Mo) 



We now determined six of the nine equations which gov- 
ern the evolution of the perturbation in the first-order 
scheme. The remaining three equations are equations 
(89), which also arise as antisymmetric part of (55), if 
we understand these as a set of nine evolution equations. 



2. Parts of the Weyl tensor m the first-order scheme 

First, we determine the electric part of the Weyl tensor 
up to the first order and find 



{i)E\ 



"E' 



P\ 



2HP\ 



AttG 



+ ^^QhP5',. (103) 



where ^iJ* is the homogeneous part. 



He^ 



-^S^ 



-{4nGpH~A)S^^ 



(104) 



which reproduces Friedmann's acceleration law (96). As 
in the Newtonian analogue we are able to rewrite the first 
part of the ADM equations in terms of the electric part 
of the Weyl tensor, i.e. 



We'' 



:k] 



(99) 
(89) 



(105) 
(106) 



On the other hand, using (89) we find the magnetic 
part of the Weyl tensor to be 



Wh^ 



1 



.ikl r> 



Ml 



(107) 



so iJ' ■ has no homogeneous part, ^H'^- = . The vanish- 
ing of its trace reproduces (89), whereas the vanishing of 
the antisymmetric part reproduces the momentum con- 
straints. 



«i/\. =0 



(i)i/. 
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(89), 
(90). 



(108) 
(109) 



At first order the relations between the parts of the 
Weyl tensor and shear tensor, respectively Ricci tensor, 
become somewhat simpler. With the first-order shear 
tensor. 



(iV% 



p\ 



\ps^, 



(110) 



we find from (71), 



(111) 



respectively ^^'^E^, — —^^'^a^, with lowered index. With 
the first-order term for the magnetic part of the Weyl 
tensor above, its relation to the time-derivative of the 
3-Ricci tensor simplifies compared with (B.2), 



2^™(,e 



mkl 



ii)H, 



-^'^R, 



nm a ^"'^ 



(112) 



3. General solution for the first-order trace part 

In this section we shall derive a general solution for 
the trace evolution equation (99). The homogeneous 
Friedmann equation (96) determines the scale factor a{t). 
With this, we separate the time and spatial derivatives 
and make the ansatz 

p-^(X,t) = °Q^(X) + q^itrQ^^iX) + q2{tfQ%{X) , 

where the time functions qi/2{t) are the two solutions of 
the linear differential equation 



2-9+(3--A)(g + g(to))=0 
a \ a 



(114) 



Note that the first part in the ansatz reflects the nonvan- 
ishing of the initial perturbation held since here we have 
to take into account initial conditions for the perturba- 
tion field as well as its first and second time derivatives. 
Solving (100) instead of (99) we have three solutions in 
a natural way since there the differential equation is in- 
homogeneous. (The result is the same as we show in the 
example in the appendix where we explicitly solve the 
inhomogeneous equation.) Note also that if we insert the 
above ansatz for P° {X, t) into the first-order Raychaud- 
huri equation (99), we get the constraint 



Ona 



\{X) = q,{t^YQ\{X)+q2{t^YQ\{X) . 



With the ansatz and its time derivatives we find 
92(to)/^",(to)-92(io)i^1(io) 



'Q\ 



2/na 



qi{to)q2{to) - qi{to)q2{h) 

qi{to)P-ito) - qi{t„)P-So) 
qi{to)q2{to) - qi{to)q2{tQ) 



and 



O/^a TDO, 



Q% = PMt^)-qi{hYQ\^q2{h) 



\2ria 



(115) 
(116) 

(117) 



Hence, the first-order "peculiar" -coframe takes the fol- 
lowing form: 

'^'''V^ = V^ + (91 (i) - qi{h)YQ% + (?2(0 - q2{t,)YQ%. 

(118) 
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where ?7° = (5° +P"(io) is the coframe at the initial time. 
FinaUy, we define the pecuhar-quantities u°- and w°'^ by 



ir\=Hi^\+u\, u^.^aP^ 



a 
They are related to the spatial functions Q\ by 

w , — 2ii + ^ m ,, = a ^ y , (119) 



91 



91 



so we are able to express the first-order coframe in terms 
of these quantities. 



IV. RELATIVISTIC ZEL'DOVICH 
APPROXIMATION (RZA) 

First, a comment concerning the use of the wording 
"Relativistic Zel'dovich Approximation" in previous pa- 
pers is in order, avoiding from the beginning of this sec- 
tion confusions that may arise during the presentation. 
Previous work, e.g. [51], suggested to use the relativis- 
tic deformation field (see (125) below) in analogy to the 
Lagrangian deformation gradient of spatial derivatives of 
the Newtonian trajectory field (21). While the density 
field is then calculated through its exact integral (35), 
evaluated for this deformation field in accord with the 
Newtonian expression (22), the spatial metric, the spa- 
tial Ricci curvature and other variables are still taken to 
be those of the strictly linearized case. In order to explain 
why our point of view will differ, we recall Zel'dovich's ex- 
trapolation idea at the basis of his approximation [76, 79- 
81] in relation to the exact foundations of this approx- 
imation in the framework of a Lagrangian perturbation 
theory [14, 15]. Zel'dovich indeed used the exact integral 
for the density field in the Lagrangian picture of fluid 
motion, well-studied in the context of an inertial contin- 
uum, and adjusted the coefficient functions in that inte- 
gral, so that its linearization on a homogeneous-isotropic 
background cosmology would reproduce the result of the 
Eulerian linear perturbation theory. He considered only 
the growing mode solution that, asymptotically, supports 
the parallelity condition (19). While in the beginning he 
combined the exact solution for an inertial continuum 
with the linearized solution of gravitational instability, it 
was later confirmed by Doroshkevich et al. [40] that this 
ansatz for the trajectory field also self-consistcntly solves 
the divergence equation for the peculiar-field strength. 
Thus, the motivation was born by appealing to the exact 
solution of the inertial continuum [82] , and Doroshkevich 
et al. (loc.cit.) added the decisive consistency test in 
the framework of self-gravitating motion. However, as 
we shall explicitly explain below, we can strictly define 
the extrapolation idea in the framework of the full set of 
Lagrangian equations for self-gravitating motion, as was 



done in [14], and we employ this definition also in the 
relativistic context. One consequence of our definition 
is that we shall functionally evaluate all field variables 
without linearizing the functional expressions. 



A. Definition: The "Relativistic Zel'dovich 
Approximation" (RZA) 

Within a fully Lagrangian perturbation approach the 
trajectory field (21) solves the full Lagrange-Newton- 
System (7,8) to first order, while it is important to em- 
phasize that this latter system exclusively contains the 
deformation gradient ff^ as the only dynamical field vari- 
able. Hence, this fact suggests that it is possible to con- 
sider the first-order solution of the LNS as an input into 
other fields that, like the density integral, are just def- 
initions and can be functionally evaluated for any tra- 
jectory field. Given this remark we aimed in the present 
work to also write down Einstein's equations for only one 
field variable, the nonintegrable deformation coefficients 
r]"'^, and consider the linearized solution as an input into 
functionals of these deformation coefficients. 

We think that this point of view best reflects 
Zel'dovich's extrapolation idea. If we would linearize all 
functionals, e.g. the metric as a quadratic form of the 
deformation field, we would just repeat the well-known 
result of the linearized Einstein equations. If we would 
keep one nonlinear expression like the density integral, 
the extrapolation idea is not consistently applied. Note 
that by taking the exact density integral we make sure 
that mass is conserved for any given perturbative solu- 
tion; but also: by taking the quadratic form of the met- 
ric we make sure that we measure distances correctly 
for any given perturbative solution (being important for 
setting up the light cone structure for a given order of 
approximation); similar remarks apply for other fields. 
We therefore propose the following definition. 

Definition: "Relativistic Zel'dovich Approximation" 

We consider the 9 functions in the coframe coefEcients 
as the only variables in the full set of ADM equations 
for the matter model "irrotational dust" within a Row- 
orthogonal foliation of spacetime. We then consider the 
general linearized solution for these coefEcients. The ap- 
proximation "RZA" consists in exactly evaluating any 
other field as a functional of the linearized solution, with- 
out performing further approximations or truncations. 

The following is a restriction that complies with the use 
of the Newtonian form of Zel'dovich's approximation. We 
restrict the general first-order solution (118) to its trace 
part and subject it to the condition 



^Q^(X) = 



(121) 



i.e. we only consider the growing mode solution. Hence, 
with (120) the Zel'dovich restriction for the peculiar- 
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fields reads 



9i- 



(122) 



or, in terms of the initial perturbation field via (116), 

(123) 



P-^(X,to) = ^P'',{X,to) 



qiito) 

With this restriction we have 

qi(to) 

and we find for the "pecuhar"-cofranie, 

-V,(X,t) = S% + P'^,iX,to)+mP''.{X,to) 
where 

_ qi{t) -qi{to) 



m 



qi{to) 



(124) 



(125) 



(126) 



Remember that the "peculiar" -coframes were defined by 

Apart from the term that arises because of the nonva- 
nishing initial perturbation, this solution is familiar from 
the section on Newtonian dynamics above. The corre- 
sponding expression for the deviation field p"j takes, be- 
cause of P",(to) = p",(io) - i?(io)p"i(to), the form: 

KZA~a^ = 5« + (1 - Hito)) p^So) + mP^ito) . 

We furthermore suggest, and we imply this in our gen- 
eral definition above, to extend the extrapolating approx- 
imation RZA to the general first-order solution, notably 
including its trace-free part where, this latter, allows to 
define a nonlinear approximation for gravitational radia- 
tion. 



B. Functional evaluation of field variables 

As a consequence of the above definition we are now 
going to evaluate a number of relevant fields. In this 
section we write 

P^,^P'^,iX,to) and P'^,^P",(X,to), (127) 

because a(t), ^(t) are the only time-dependent functions. 
As said in the above definition of the RZA, we insert 
the coframe ^^'^■q°-^ = air}"-^ + ^P"), where 



f,\^5^, + P'^,{XM) 



(128) 



is the coframe at some initial time to, into the exact def- 
initions and equations, as given in section III. The sym- 
metry condition (62), for example, is extrapolated to 



{2Hi + S,){Pu.+5abP\P\)^Q 



(129) 



The metric for the coframe ^^^■q'^. — 0(77° + ^-P" ) 
takes the following quadratic form: 

"""% = a\, [fi'^^f,'^ + 2^ri^P% + eP^.^P)] ■ (130) 
Furthermore, we define: 



•'' j — c^ahc^ V jV kV I ) 

ikl pa 'b 'c 
X^abc^ ^ jV kV I 



(1)J\- 



ikl °a r)b °c 



(3)j. 



, ild'a po o. 



I I 



ikl pa pb »c 
^abc^ ^ j^ k'l I 



ikl pa pb pc 



ikl °a pb 



6 



iKl'a pO p 



I ' 



With the help of these definitions the transformation 
determinant "^*J reads: 



^^\J^a\j^+£,J,+eJ2+eJz), 



(131) 



where J„ = ^"1/*^. In the homogeneous case, i.e. for van- 
ishing initial perturbations, we find ^Jq = 1 whereas ^Ji- 
^J^ reduce to the three scalar invariants for P". Hence, 
the nonlinearly evolved density can be calculated through 
its exact integral, 



-.^=_^ = ^^(X) ^ 



""V QHito 



J 



(132) 



where '^^V ee det(fy''J = Jq + ^-^1 + ^^-^2 + C'Vg is eval- 
uated for the "peculiar" -coframe field (125) and J = 
J(A', to) is the transformation determinant at some ini- 
tial time. 

Then the orthonormal vector basis e^' in the RZA- 
picture is given by 

so the orthonormality relation for 77°^ and e^' in the RZA 
becomes (here and in the following the summation always 
runs from n = . . . 3 ) : 



RZA I RZA a 

^a '/ i 






with trace ^'^^ej^ ^'^^if^ ~ 3, as expected. 

To be able to write the RZA expansion tensor in a 
similarly short form we define the second set of auxiliary 
quantities, i.e. 



(0)p 



Wn = 



{2)p 
{3)p 



1 

6' 
1 

2 
1 

2 

ikl pa pb pc 
'a^abc^ ^ j^ fc-' I ■ 



Akl jLa °b °c 



ikl pa 06 oc 
^abc^ ^ jV kV I J 



ikl pa pb »c 
^abc^ ^ j^ kn I ; 
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Of course '^°¥j = (°l/^^ as well as <-^¥^ = '•^^J'j and their 
traces are identical to those of the J'^j's. 

Hence, the RZA expansion tensor takes the form: 



^Q^ 



3H 






(133) 



We evaluate the mixed components of the tensor, since 
in this form the corresponding Newtonian expressions arc 
easily recovered. The expansion scalar, for example, is 



^e = 3iJ + T- 



i ^71+2^^2 + 3^^3 



^ jo+^j,+ej2+ej3 



(134) 



Now wc give the expressions for the parts of the Weyl 
tensor in the RZA. The electric part reads: 



3v i^zAj ; J 



(135) 



The relation to the mixed shear tensor components reads: 



A^^ 



RZA^i RZA^A; 



2^ En<"4 

3^ E„e4 



RZA „k RZA ^/ r2 



(136) 



(This should be compared with the somewhat simpler 
form (111) in the first-order scheme, and with the general 
relation (71).) 

The magnetic part takes the following form: 



"-W 



_^^ikl , fejil' 



Kb{np\+^np' 



k)\\i 



E„e^„ 



(137) 

where || here denotes the covariant derivative with re- 
spect to the RZA-metric (130). 

Finally, we express the Riemann and Ricci curvature 
tensors functionally in terms of the RZA-deformation. 
To keep the equations short, we define (analogous to the 
functionals ^°l7' ■ — ^'^-!/' above) the quantities 



(")ff 



jkl — p^afcc^ V j\\k\\lV mn n i 

1 



^ jkl = T-^abc^ ^ j\\k\\lV mV 



(2)^» 



jkl 



imn'a r,b »c 



b 'c 



- imn pa pD =. 

l^abc'^ ^ ]\\k\\l^ mV 



I imn 'a pb pc 

+ p^afcc^ V j\\k\\l-^ m^ n i 

^ imn pa pb pc 

T-^abc^ ^ j\\k\\l-^ m-^ n ■ 

The curvature tensor then is 



(3)^. 



RZApZ o 

^ jkl — -=> 



En S 'Jn 



(138) 



and the Ricci tensor takes the form 



X 



E„r('"^^' 



ikj 



{n)Rl 



ijk) 



En S 'Jn 



(139) 



Of course one can express the curvature quantities solely 
in terms of the RZA-deformation, since we used the co- 
variant derivative with respect to the RZA-metric (130), 
but the so-found equations are long and hard to read. 



V. DISCUSSION AND CONCLUSIONS 

In this section we put some aspects into perspective 
that were discussed in the context of the Newtonian 
"Zel'dovich approximation" , and we summarize the main 
points of this paper. 



A. Possible tests of the extrapolation 

Applications of the presented approximate theory leave 
in suspense the unknown quality of the extrapolation 
done. Of course, we may blindly accept the resulting 
approximate solutions driven by the belief that the cor- 
responding Newtonian model is extremely successful in 
comparison with N-body simulations of the full prob- 
lem. However, it is in order to point out that a num- 
ber of self-consistency tests are possible and should be 
done. While this scheme can predict and describe ef- 
fects beyond the known highly symmetric solutions of 
general relativity, it is necessary to conduct additional 
tests. Such tests are best performed in the context of 
the envisaged applications, and we shall come back to 
them in forthcoming papers. Recall that Doroshkevich 
et al. [40] provided such a test for the Newtonian form of 
Zel'dovich's approximation by considering quantitatively 
the difference between the density calculated from the 
exact integral of the continuity equation and the density 
calculated from the field equation, i.e. from the diver- 
gence of the peculiar-gravitational field strength. While 
at first order both expressions agree by construction, the 
error was reported to be of second and higher order [40] . 
While qualitatively this is obvious, the error was calcu- 
lated and used to estimate the quantitative validity of the 
approximation in certain regimes. In this spirit we can 
also compare resulting nonlinear expressions and conduct 
consistency tests. As an example we note that e.g. the 
scalar curvature can be determined from the Hamilton 
constraint (29) through the kinematical invariants and 
the density in the RZA, and alternatively from the RZA 
metric by explicitly calculating the trace of its Ricci ten- 
sor (53). The resulting expressions agree to first order 
and the error is of higher order and may be quantita- 
tively controlled in the context of a given application. 
We shall come back to these different curvature expres- 
sions in forthcoming papers. 
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B. Singularities 

Continua made of dust are bound to develop singular- 
ities in the course of evolution, resulting in caustics, i.e. 
loci of formally infinite density. This is a consequence 
of neglecting physical effects like velocity dispersion or 
vorticity that could regularize singularities (for the New- 
tonian theory see [28] ) . In general relativity this situation 
corresponds to the intersection of world-lines and, thus, 
to the failure of defining a congruence of world-lines, to- 
gether with the possibility of simultaneously developing 
singularities in geometrical fields. 

The appropriate mathematical framework in which 
caustics can be described and classified is catastro- 
phe theory [77], further developed in the framework of 
the Lagrange-singularity theory especially by Vladimir 
Arnol'd and collaborators. Caustics are defined as images 
of a critical set on a Lagrangian submanifold under a pro- 
jection map. The stable singularities of such Lagrangian 
mappings can be classified into a finite set of topologi- 
cal structures (germs), their number depending on the 
number of dimensions of the Lagrangian submanifold. A 
classification into a finite number of elements is only pos- 
sible for manifolds with dimension < 5. Alternatively, 
the Lagrange-singularities of a family of world-lines of 
fiuid elements can be described in terms of Legendre- 
singularities of the wavefronts (if they exist) as the dual 
description of the continuum's evolution (c/. the small 
selection of papers by Arnol'd [3], [4], [5], the books [7] 
and the collection of papers [8]). 

The singularities developing in an irrotational New- 
tonian continuum moving under inertia have been ana- 
lyzed in full detail with nice hand drawings by Vladimir 
Arnol'd [3] , and a formal relation of this mapping to the 
(Newtonian) "Zel'dovich approximation" has been estab- 
lished and analyzed in detail for caustics in two spatial 
dimensions [9]. In the pancake picture developed at the 
time one considers images of singularities as the local 
building blocks of the large-scale structure in the Uni- 
verse. The geometry of these structures may differ for 
different realizations of the model, but the local mor- 
phology of structures (so-called unfoldings around sin- 
gularities and their evolution (so-called metamorphoses) 
is completely made up from 12 (topologically classified) 
elements in a four-dimensional continuum [3]. 

Specifically, in Newtonian theory, we may define a flow 
field X — f (X, t) by a time-dependent diffeomorphism 
that sends initial (Lagrangian) positions X of fluid ele- 
ments to their Eulerian positions x (embedded into Eu- 
clidian space) at time t. In order to apply the clas- 
sification scheme of the Lagrange singularity theory to 
the motion of a Newtonian continuum, the key-property 
that has to be required for f is that it can be written 
as a family of gradient mappings, which form an impor- 
tant class of Lagrangian mappings. This is for example 
easily possible, if the flow is irrotational with respect to 
the Lagrangian frame, i.e., there exists a potential ^ for 
which f —: Vx^, where Vx denotes derivative with re- 



spect to Lagrangian coordinates. We may then define 
the one-parameter family (parametrized by the time t) 
of Lagrangian mappings nt : 



Tl-t 



Xh^x = Vx*(X;f) 



(140) 



where in this case the set {X,x = f(X, t)} forms a La- 
grangian submanifold of MP = {X,x}. Then, for each 
fixed t, the (non-degenerate closed) two-form J2i ^^^ ^ 
dX^ vanishes on the Lagrangian submanifold (X* are 
local coordinates on this manifold). Note that the re- 
quirement of irrotationality of the flow field with re- 
spect to Lagrangian coordinates is in general much more 
restrictive than the requirement of irrotationality with 
respect to Eulerian coordinates. For rotational flows, 
the Lagrange-singularity theory has been extended by 
Bruce and collaborators [13], specifically investigated for 
solenoidal velocity fields of an inertial continuum. They 
found that the most ubiquitous pancafce-singularities {A3 
in Arnol'd's classification) remain stable. However, the 
singularities associated with umbilics in the potential 
case have to be removed from the list of generic (struc- 
turally stable) patterns in the vortical case. 

The possibility of writing the flow fleld f in the form of 
a gradient mapping can be demonstrated for some impor- 
tant subclasses of irrotational Lagrangian perturbation 
solutions: the flow fields develop Lagrange-singularites in 
the case of first-order solutions [15], and for a large class 
of second-order solutions [29]. The third- and fourth- 
order contributions destroy this possibility even for re- 
stricted classes of initial conditions [17, 70]. 

In the case of self-gravitating continua, a proof of the 
Lagrangian property meets the problem that the veloc- 
ity field and the acceleration field may become multi- 
valued simultaneously ([6] §6, footnote 1). Nonperturba- 
tively the situation is worse, since the gravitational field 
strength does not remain finite at caustics, as would be 
suggested by Lagrangian perturbation solutions to any 
order [22, 23]. Indeed, a consequence of the action of 
self-gravity is the development of a hierarchy of nested 
caustics, which originate simultaneously, with an ever in- 
creasing number of streams (see [29] for an example in 
second-order perturbation solutions). 

In general relativity the problems described above are 
more involved. One open problem is the dynamical jus- 
tification of the continuation of solutions across caus- 
tics, where the multi-stream flow is required to satisfy 
the field equations. A mathematically well-defined con- 
cept of bifurcating dust appears to be a difficult prob- 
lem; it is evident that also the shapes of caustic sur- 
faces may not necessarily permit a morphological clas- 
sification in the framework of the Lagrange-singularity 
theory. Clarke and O'Donnel [37] succeeded in showing 
the self-consistency of an extension of spacetime through 
a dust caustic. Other authors address the problem in the 
spherically symmetric case [50]. In cosmology the sin- 
gularity problem is mostly addressed for the asymptotic 
past of solutions; for a study of the asymptotic structure 
of cosmological singularities see, e.g., [2, 58]. Of course 
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there is a substantial literature related to singularities in 
general relativity (see, e.g., [36, 55, 66], and references 
therein, as well as the book by Clarke [35]. 

The model proposed in this paper contributes to this 
discussion. We can explicitly study those fields that 
would feature singularities in the relativistic Zel'dovich 
approximation. Since caustics correspond, in the New- 
tonian model, to degeneracies of the Jacobian, J — 
det(/*|^,) = 0, we have to look at degeneracies of the 
corresponding relativistic field (denoted by the same let- 
ter in the present paper), J = det{r]\) = 0, in the local 
exact basis dX^. From this it is already evident that 
the nonintegrability of Cartan's deformation field - in 
general - destroys the possibility of defining a family of 
gradient mappings, and a morphological classification in 
the classical framework is not straightforward. With the 
reasonable assumption that local deformations remain fi- 
nite we conclude that the metric, as a quadratic form of 
the deformations, remains finite at caustics. However, 
as the explicit formulae in Appendix C show, almost 
all relevant fields will degenerate simultaneously, since 
they are weighted by "^^^J (while the remaining terms 
remain finite for finite deformations). Note here that 
the verticality (with respect to Eulerian coordinates) of 
the Newtonian velocity gradient, which is indicative for 
the occurrence of a singularity, corresponds to the rel- 
ativistic field Q^ that becomes the mixed-index object 
0'i = e^rf^Q^ in the exact basis, c/. Eq. (C.16). The 
explicit expressions in the RZA model imply that spa- 
tial but also spacetime curvature terms become singular 
at caustics. Due to the form of the general expressions, 
however, we expect these degeneracies to appear also in 
general situations. 



C. Concluding Remarks 

Following the systematic derivation of Zel'dovich's ap- 
proximation in the Lagrangian framework of the New- 
tonian equations we have formulated the Einstein equa- 
tions for the matter model 'irrotational dust' in terms 
of a single dynamical variable. The nine coframe coeffi- 
cients of Cartan's deformation one-forms replace the in- 
tegrable Newtonian deformation gradient deriving from 
three vector components. We discussed the resulting sys- 
tem using different representations, and we derived the 
general first-order solution for the coframe functions (the 
deformation field) . We then gave a definition of a nonper- 
turbative approximation scheme that proposes to func- 
tionally evaluate dynamical field variables in terms of 
the perturbed deformation field. The success of the cor- 
responding Newtonian approximation gives substantial 
motivation for this generalization to relativistic cosmol- 
ogy. 

There are a number of aspects that we consider use- 
ful. First, using the proposed equations we can eas- 
ily translate "Newtonian knowledge" to the relativistic 
stage. This is especially due to the formally close corre- 



spondence of the 'electric part of the LES' to the LNS of 
Newtonian theory. Second, we can employ explicit forms 
of nonperturbative expressions for field variables that just 
depend on initial data and known time-dependent func- 
tions. We so are able to attack highly nonlinear problems 
in relativistic cosmology. For example, the approximate 
quadratic form of the metric can be used to realistically 
evaluate distance expressions in inhomogeneous structure 
distributions; the explicit structure of the light cone al- 
lows the study of the influence of generic inhomogeneities. 
The same is true for the Ricci curvature and the Weyl 
curvature with its parts. As a consequence this approxi- 
mation allows to investigate many problems beyond the 
Newtonian approximation such as light propagation and 
gravitational radiation. 

In forthcoming work we extend this model by employ- 
ing exact integral properties of Einstein's equations [19]. 
The combination of a generic model for inhomogeneous 
deformations with exact integral properties has led to a 
powerful model in the Newtonian approximation [32, 52], 
and so we shall investigate the corresponding relativis- 
tic problem in order to quantify the influence of inho- 
mogeneities on average properties of the Universe (the 
'backreaction problem' [24, 25, 45]). Furthermore, we 
shall give the general perturbation and solution schemes 
for the Lagrangc-Einstcin-System including examples of 
higher-order Lagrangian perturbation solutions. Non- 
perturbative investigations of light propagation in inho- 
mogeneous models as well as nonlinear aspects of gravita- 
tional radiation are also envisaged in future applications. 

A further, more challenging but possible, application of 
the presented formalism would employ a self-consistently 
evolving background rather than a flxed FLRW back- 
ground as in our examples. A background including 
backreaction effects could be determined by the exact 
average properties of an inhomogeneous universe model 
[71]. Including pressure by translating Newtonian results 
is possible [1, 28]. Furthermore, a nonvanishing shift 
vector field together with a non-constant lapse function 
could be included and herewith the Lagrangian condition 
extended, all providing more general frameworks in the 
spirit of this work. 
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APPENDIX A: Formulation using Differential 
Forms 



Newtonian equations 

An alternative, somewhat simpler form of the 
Lagrange-Newton-System of equations can be repro- 
duced by introducing differential forms. Using the spatial 
exterior derivative operator d, acting on functions and 
forms, and the spatial exterior product, the LNS takes 
the form [42]: 



S^jdfAdP = 0: 



(A.l) 



1 



e,,fcdr A df A 01/*^ - (A-47rGe)dV, (A.2) 



where the density is given by the integral g = qJ~^, 
and d'^/ = Jd^X. Defining the three differential one- 
forms of the tidal tensor £* = fV.dA'^, the LNS assumes 
instead the compact form: 

5^j£' Adf = 0; 



where 



e„-fe5^ A dp A dp = , 



5' = df - 3 (A - AttGq) dp 



(A.3) 
(A.4) 

(A.5) 



Einstein equations 

We introduced Cartan's coframcs 77° — ?/°jdX*, 
one-forms that define a noncoordinate basis of three- 
dimensional space. We may call them spatial deformation 
one-forms. The metric form is the canonical quadratic 
form 



g = Sab-n" ® rf 



(A.6) 



Noncoordinate indices are raised and lowered by (5^j,. 

The expansion one-form is then defined by the parallel 
transport equation, i.e. 



r(\dX^ 



(A.7) 



Spatial derivatives that take into account the varying ge- 
ometry are evaluated, for a symmetric connection, by 
Cartan's connection one-forms. 



b = 7 c6^ 



(A.8) 



together with the curvature two-form 



(A.9) 



These objects are defined by the (spatial) Cartan struc- 
ture equations: 



drf + u;% hrf ^T" 



dw". + u;°„ A uj" 



n" 



b ' 



(A.IO) 
(A.ll) 
(A.12) 



together with the integrability conditions ddrf = and 
ddu;"j = 0, where T° = \T°-bjf' A rf is the torsion two- 
form. We set the torsion to zero in this work. Hence, 
both Lx>°'^ and il";, can, in principle, be expressed solely 
in terms of the deformation one-forms. 

We may define the total exterior derivative of a tensor- 
valued differential form $ by 



(D$)% = d*-^, + ^\ A -i\ - L,\ A ^\ 



(A.13) 



which is the natural generalization of the covariant 
derivative when working in a noncoordinate basis. With 
this definition we have: 



Dj7° 



DT° = n^f, A rf and D^^b 







(A.14) 
where the two last equations represent the two Bianchi 
identities. 

The invariant volume element in the noncoordinate ba- 
sis is 

J7I A 77^ A r/3 == JdX^ A dX^ A dX^ = Jd^X . (A.15) 

The Lagrange-Einstein-System takes the form 

<5a6J7° A r?^ = , (A. 16) 

(A J - 47rG^) d^AT , (A.17) 

, (A.18) 






and 
1 



-ebcdV A r/'- A J7" 



ebcdil" Ml' hv'' (A.19) 

[{AttGq + A J) (Jf," - JT^b"] d^^ , 



where the 3-Ricci tensor can be expressed via the curva- 
ture two-form by 

W^ri'^ Ari'' A rf' = S'^'^Q,"^ A 77== - S'^'^Q," ^ A rf . (A.20) 



Let $ = <i>ij...i^da;*i A... A da;*'' be an r-form in a three- 
dimensional manifold and g = dct{gij) the determinant 
of the metric. Then the duality operator, or Hodge star 
operator, is defined by 



* $ 



(3-0 



■$,: 



r+l---«3 



dx' 



A dx^^ . 
(A.21) 
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In particular, *1 is the invariant volume element because 
of 



*1 = 



9 



0| '^»1»2*3'- 



A dx'^ A dx'^ 



(A.22) 



With the help of (A. 21) we can also introduce one- 
forms of the electric (tidal) and magnetic parts of the 
Weyl two-form, which are 



which together with (33) allows us to relate the pro- 
jected magnetic part of the Weyl tensor back to the time- 
derivative of the 3-Ricci tensor. 



R, 



-25. 



m(i-'" ^J)fc||; +0y \\k-^\\i\\j 



(B.2) 



To rewrite (34) in terms of the Weyl tensor we now de- 
rive a full set of evolution equations for the electric and 
magnetic parts. 



E" 



»7 



(A - At:Gq) r," 



H« = *(<5,,Dr,^A77")r,'' 



(A.23) 
(A.24) 



The differential form counterparts to the equations (74) 
- (77) then read: 



eab^E" A 77^ A T7= = : 

5„bH° A r,'' = : 

e„bcH" A ry^ A r?^ = . 



(A.25) 
(A.26) 
(A.27) 
(A.28) 



We see that the 'electric part' of Einstein's equations, 
(A. 16, A. 17), or (A.25, A.26), respectively, immediately 
reduce to those of Newtonian gravitation (up to the sign 
convention) for exact differential forms as geometrical 
limits of the coframes, i.e., 






r\kdX' 



(A.29) 



Applying this geometrical limit to the magnetic one- 
form (A.24), we immediately see that it vanishes in the 
Newtonian limit, since the coefficients of the expansion 
one-forms 0° = r/" reduce to the velocity gradient 



NQa _^ p dX'' ^ v^'i^dx'' with partial derivatives. 



^H" ^ S,f'""^v''^Ax'' = 0. 



(A.30) 



where the Newtonian limit implies the existence of global 
Eulerian coordinates x''. (For a useful reference on Car- 
tan's formalism see [65].) 



APPENDIX B: Other formulations involving 
electric and magnetic parts of the Weyl tensor 

We here provide other formulations in terms of the 
Weyl tensor parts that are helpful especially in future 
work. 

Note that throughout this appendix we adopt the in- 
variant volume element e*'^' rather than the antisymmet- 
ric tensor e**^', see (45). 

An interesting identity for the magnetic part is a rela- 
tion of its covariant curl to covariant spatial derivatives 
of the expansion tensor. 



9^: 



'H. 



3k\\l 



= e 



\\k 

ij II fc 



-e' 



3M\k 



(B.l) 



Maxwell— like equations for the projected parts of 
the Weyl tensor 

The electric and magnetic parts (68) of the Weyl tensor 
obey a set of evolution equations similar to the Maxwell 
equations, see for example [46, 78], [11, 47, 53, 59]. The 
Weyl tensor plays the role of the electrodynamical field 
tensor. Starting from the 4-Bianchi identities in the case 
of irrotational dust, u^ = (1,0,0,0) and w^.^ ~ Q^^ 
(where a semicolon denotes 4— covariant derivative). 



C 



UKX.fl 



8ttG 



(qu^'^Ix); 



1 



:0;lK9,y\] 



(B.3) 



we replace the Weyl tensor by its parts (68), 



Cf^^KX = {9,„.ap9KXjS - ^^.i.ap^KX'is) u°'u'^E^^ (B.4) 

where g^^„^ = 9f,a9^fi - g^f)9ua and the te nsor of the 
invariant volume element s -^ ~ v ^'"^''s^ui/kA ^^ '^^~ 
fined analogous to (45). We transform to Lagrangian 
coordinates as before and find, with y— ^^5 = y^ = J 



and £ 



fil^fiX 



^ukX ^°^ ^^^ spatial invariant volume ele- 
ment, the following equations for the time-derivatives of 
the projected parts of the Weyl tensor or their covariant 
curls, respectively: 



E' 



W 



2QE'' 



QjE\ 



Q'',E\,5\ - e'''^ H 



i^ fe" ] 



Mk 



= -4-G-^(e^,. 



ie.^ 



(B.5) 



2ew 

A-kG J 



e'^W, Q\H\5^^ 



- £*'='K 



Mk 



— J9\ke 5,.. 



(B.6) 



Recalling that taking the time-derivative does not com- 
mute with the raising and lowering of indices, i.e.. 



E., - i9^kE',) - 9^kE', + 2e,,E'^ 



(B.7) 



(analogous for the magnetic part), and splitting the ex- 
pansion tensor into its kinematical parts (for vanishing 
vorticity). 



6' 



(B.8) 
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wc get the well-known forms of (B.5), (B.6) as given in 
the literature, e.g. [11], [47]. For the covariant diver- 
gences of the Weyl tensor parts wc find: 



k „ c^kmnr\ rrl 

i\\k ^ 9ik^ '^ml^ n 



SttGJ 



3 J^l* ' 



^'||fc+5.fe^'™"0w^'„ = O 



(B.9) 
(B.IO) 



The antisymmetric parts of Equations (B.5) and (B.6) 
are equivalent to Equations (B.9) and (B.IO), so we just 
take their symmetric parts (with lowered indices) into 
account, i.e., 



E, 



•y I 2e£;,y - SQ^iE j^ -& lE ^.g^j - g^^iS^ ^m\k 



.J. 
J 



1 



-47rG-^(ey--e5,, 



(B.ll) 



Hij + 20-ff.y - 36j,(jiJ ^-^ -Q iH f.gij + .g„(je"' £'j)i||fc 
= 0. (B.12) 

Note that for the right-hand-side of (B.12) we applied 



AttG J 



%kl 



3 J 



T^lfeE 5r„(»5,)i - 



(B.13) 



which is due to the antisymmetry of the Levi-Civita- 
tensor density. These equations are automatically satis- 
fied, if we insert the expressions for the parts of the Weyl 
tensor we found in the previous subsection and apply the 
ADM equations. Equation (B.ll) is the equivalent to 
(34) we looked for. 



APPENDIX C: EXAMPLE FOR THE RZA 

In this appendix we shall derive a general solution for 
the trace evolution equation (100) for a flat FLRW back- 
ground universe. (This way we also clarify the way in 
which the trace equations (99) and (100) are equivalent.) 
The homogeneous Friedmann equation that we restrict to 
the Einstein-de Sitter background, determines the scale 
factor 



The time function q{t) then has to obey the inhomogc- 
neous differential equation 



q- 



2 . 1 /.., N 2 ., 
79= —(9^0) + — 9 t 



(C.3) 



Since (C.3) is linear and of second order, the general so- 
lution can be written as superposition of three solutions, 
so a = 1,2, p. (Here p denotes the particular solution of 
the inhomogeneous equations, whereas a = 1,2 label the 
homogeneous solutions.) With the ansatz 



we find 



Hence, 





9V2W = (J ' '^^(^)-U) 


id 


2 
ni = 0, 712 = -1 , -P= 3 • 


■^5 


^■Q". + (i)"V. + (i;;)''<3V 



(C.4) 



(C.5) 



Using this we express the initial perturbation fields in 
terms of the initial conditions P° (io), P°'iito) and P° (io), 



^-P-,{to)4~2P-So)to + P''dto). 



2/na 



pn°-. 



^p-^{to)tl + ^P'^,{to)to. 



Thus, the general solution to (C.3) is, written in terms 
of the coframc for the deviation field: 



«77", = 5", + P'^,(io) 
1 / i \-i 



(C.6) 



lOVto 



10\toJ 



2/t\-i 3/<\t 






«(^) - a 



(C.l) 



The RZA for a flat FLRW background with A = 



After solving the equations we restrict the solution to the 
Zel'dovich approximation and give some of the quantities 
discussed above for this particular choice of model. 



As in section IV we write 



P'', = P-,{X,to) and P'', = P'',{X,to) 



General solution for the first— order trace part 

For the inhomogeneous equation (100) we make the 
ansatz 



p^,ix,t) = J2i<^it)"QM^)' 



(C.2) 



because a{t), ^(t) are the only time-dependent functions. 
First, we determine Zel'dovich's restriction in both 
forms (122) and (123) for the flat FLRW background with 
vanishing cosmological constant (Einstein-de Sitter back- 
ground). The growing mode in (C.5) is q{t) = (i/ig)^^'^, 
so with the scale factor (C.l) we flnd 



<(to)=<(to)to, 



(C.7) 
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and, in terms of the initial perturbation fields, 



(C.8) 



We subject the general first-order solution (C.6) to 
Zel'dovich's restriction and find, with 



tensor take the following forms: 



RZAq^I^(1\-^ EnC-'Jn 



t Vt, 



RZAJ 



and 



(CIS) 



«0 ^ I 



b--^ 



for the "peculiar" -coframe: 



^7?", = 5'^, + P'^+2 






p\t^ 



(C.9) 



(CIO) 



Apart from the term that arises because of the nonvan- 
ishing initial perturbation, this solution is familiar from 
the section on Newtonian dynamics above. In the follow- 
ing we shall use ^ and the ^"l/* ■ of the RZA section to 
keep the equations short. 

Since ff\ = (5"^ + P° , we write 

Hence, the metric coefficients take the following form: 



9'. 



y 






-]"K^i^i\n', + -^in^fi"",) + ep\p\) 



and the transformation determinant ''^*J reads: 

3 



(C.12) 



J = (-) "^V , with -V ^ Y. ^"-^^ ' (C.13) 



n=0 



where J„ = ^"^J\ and J = det(77"j) is evaluated for the 
"peculiar" -coframe field (C.IO). 

As in the general case, the nonlinearly evolved density 
can be calculated through its exact integral. 



-A, = _^ ^ Wg) .(^) J 



"^V QH{to) 



J 



(C.14) 



For the flat FLRW background the RZA expansion 
scalar and the mixed components of the RZA expansion 



ie» 



t 



RZAJ V f 



AJ 



(C.16) 

(Here and in the following the summation is over n ~ 
0...3.) 

The Riemann curvature tensor (138) then becomes: 



RZAp« Q 

^ jkl ~ •J ■ 



ErC"^i?V,-("^i?^- 



Ik) 



KZAJ 



and the Ricci tensor (139) takes the form: 



RZAp o 



RZAJ 



(C.17) 



(CIS) 



Finally, we give the expressions for the parts of the 
Weyl tensor in the RZA. First, with 



H ^ 



3i 



e + 2ffe = 



1 /t\-l 



to ^to 



(C.19) 



the electric part reads: 



fl,ZA£,._A.?^"^"^"^^^ 47rG gnit) pJ 



3t2 



AJ 



.S' 



id 

and the magnetic part: 



to Vio 



3 QH{to)''^M ^ 
^ (C.20) 



ikl 



RZAjji ^ _l!jl 



AJ 



Pkm+^abihP'k + ihP' 



J kj\li 



RZAJ 



(C.21) 

where || is the covariant derivative with respect to the 
RZA-metric (C.12). 
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